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ABSTRACT 

We an alyze the cosmological simulations performed in the recent work of iGreif et al.l 
(|2012t) . which followed the early growth and merger history of Pop III stars while 
resolving scales as small as 0.05 R©. This is the first set of cosmological simulations to 
self-consistently resolve the rotation and internal structure of Pop III protostars. We 
find that Pop III stars form under significant rotational support which is maintained 
for the duration of the simulations. The protostellar surfaces spin from ^50% to nearly 
100% of Keplcrian rotational velocity. These rotation rates persist after experiencing 
multiple stellar merger events. In the brief time period simulated (~ 10 yr), the proto- 
stars show little indication of convective instability, and their properties furthermore 
show little correlation with the properties of their host minihaloes. If Pop III proto- 
stars within this range of environments generally form with high degrees of rotational 
support, and if this rotational support is maintained for a sufficient amount of time, 
this has a number of crucial implications for Pop III evolution and nucleosynthesis, 
as well as the possibility for Pop III pair-instability supernovae, and the question of 
whether the first stars produced gamma-ray bursts. 
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1 INTRODUCTION 

The formation of the first stars marked a crucial transition 
in th e evolution of the early Universe (|Bromm fc Larson] 
12004 lLoebl l20ld ). These stars were metal-free, and are 
thus also known as Population III (Pop III). They are 
believed to have formed at z > 20 with in dark mat- 
ter minihaloes of mass ~ 10 6 M (e.g . Haiman et al.l 
1 19961 : iTegmark et~aH Il997l ; lYoshida et aLTl2003l ). Pop III 
stars that were sufficiently massive then formed the first 
Hn regions, beginning, the process of reio nization (e.g. 



cgirn 

Kitayama et al.l 12004 iSokasian et all 1200-4 IWhalen et al.l 
20041 : [Alvarez et al.ll2006l : I Johnson et al.ll2007l ). Stars within 



certain mass ranges ended their lives as supernovae, thereby 
contributing to the initial enrichment of the intergalactic 
medium (IGM) with heavy elements dMadau et all l200ll; 
Mori et al.1 120021: iBromm et alj 120031: IWada fc Venkatesan 



20031 : iNorman et al.ll2004lTornatore et al.ll2007l : IGreif et a! 



l2007l . l2010l:IWise fc Abel 20081: iMaio et alj|201ll : recently re- 
viewed in lKarlsson et al. |2012| ). In particular, non-rotating 
primordial stars with mass 140 Mq < M* < 260 Mq are 
believed to have exploded as pair-instability supernovae 
(PISNe: IfTeger fc Wooslevl |2002| ). releasing the entirety of 
their metal content into the IGM and surrounding haloes, 
while stars within the range 15 Mq < M* < 40 Mq ended 
their lives as core-collapse SNe. On the other hand, non- 
rotating Pop III stars with main sequence masses in the 
range 40 M Q < M, < 140 M or M„ > 260 M are ex- 
pected to collapse directly into black holes, thus contribut- 
ing no metals to their surroundings. 



The typical mass of Pop III stars is therefore a crucial 
factor in determining their role in the evolution of the early 
Universe. However, the rotation of a star ca n also strongly 
influence stellar ev olution and death, (e.g. IMaeder 1987 : 
see also reviews by IMaeder fc Mevnetl |2000| . 120121 : lLangerl 
20121). with stron ger effects for lower-metallicity stars (e.g 



Brott et al.ll201ll ). For instance, instead of a normal core- 
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collapse SN, rapidly rotating stars may end their lives as 
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extremely energetic hypernovae (e.g. lNomoto et al"1l2003l ). A 
further consequence may be that rapid rotation lowers the 
minimum mass at which Pop III stars can beco me PISNe 
(|Chatzopoulos fc Wheelerll2012l ; lYoon et alj|2012j ), possibly 
to masses as low as 65 Mq. This occurs through the pro- 
cess of rotationally induced mixing, which in turn leads to 
chemically homogeneous evolution (CHE) and a larger final 
oxygen core mass as compared with a non-rotating star of 
the same initial mass. 

Theoretical studies of stellar evolution furthermore find 
that, compared to non-rotating stars of the same mass, the 
typical effective temperatures and luminosities of rotating 
Pop III stars are larger, and for an extended ran ge of masses 
CHE shifts stellar evolutionary tracks blueward (|Yoon et al.l 
l2012h . This depends sensitively, however, on how angular 
momentum redistribution by magnetic fields is modeled. 
Earlier work which d oes not implement the Spruit-Taylor 
dynamo (|Spruitll2002h , for instance, find that Pop III stars 
will not undergo CHE and will in fact end their evolution in a 
redder part of the Hertzsprung- Russell diagram (Ekstrom et 
al. 2008a). The metal production within rotating primordial 
stars is generally higher, and rapidly-rotating models can 
lead to 14 N yields that are several orders of magnitude larger 
than for corresponding non -rotating or slowly -rotating mod- 
els (Ekstrom et al. 2008a; Yoon et al.|[2012h . However, we 
also note that iHeger fc Wooslevl (|2010h find significant 14 N 



yields even in their non-rotating models. 

Rotation may additionally determine a Pop III star's 
potential for producing a gamma-ray burst (GRB), partic- 
ularly given the connection between lo ng-duration GRBs 
and the death of massive stars (e.g. IWooslev fc BloorrJ 
2006). This would provide a promising method of di- 
rectly probing Pop III stars in their final stages, pro- 
vided that Pop III GRBs occurred with sufficient fre- 



quencv ilBromm fc Loebl |2002|. 120061; iGou et all |2004 



iBelczvnski et all 120071 ; iNaoz fc Bromberd 120071 ). The col- 
lapsar model of GRB formation requires the presence of 
sufficient angular momentum in the progenitor for an ac- 
cretion torus t o form around the remnan t black hole (e.g. 
I Wooslevl [199^ ; [Lee fc Ramirez- Ruij|2006l ). Furthermore, to 
enable the escape of the accompanying relativistic jet from 
the s tar, the progenitor must also lose its hydrogen e nvelope 
(e.g. IZhang et all 12004 but see ISuwa fc Iokall201lh . How- 
ever, this latter requirement poses a difficulty for single-star 
progenitors, because removal of the hydrogen envelope will 
also lead to a decrease of angular momentum in the core (e.g. 
ISpruitll2002l ; IHeger et al.ll2005l ; IPetrovic et ai1l2005l ) . On the 
other hand, both conditions for collapsar GRB formation 
may be met in the case of a close binary system undergoin g 
Roche lobe overflow (e.g. lLee et aI1l2002l ; lizard et alj|2004l . 

Another route to GRB formation arises for Pop III stars 
with sufficient spin, due to the effects of rapid rotation on a 
star's nucleosynthe sis as well as its evolution off of the main 
sequence (MS; e.g. lYoon fc Langerll2005l; IWooslev fc Hegerl 



l200d:lYoon et alj2012h . For i nstance , lYoon fc Langerl (|2005l ) 
and lWooslev fc Hegerl (|2006f ) find that low-metallicity mas- 
sive stars (> 20 Mq) with rotation rates above ~ 40-50% 
of their breakup speeds may undergo rotationally induced 
mixing and CHE, allowing the star to bypass the red giant 
phase and become a Wolf-Rayet (WR) star. This may al- 
low the star to retain enough angular momentum to become 
a GRB, particularly if the star also undergoes little mass 



loss through stell ar winds, as is e x pected for low and zero- 
metallicity stars (|Kudritzkil 120021 1. lYoon et all (|2012h find 
that this route to GRB production through CHE applies to 
metal-free stars as well. 

A full understanding of the nature of Pop III stars 
and their impact on the IGM thus requires knowledge of 
their typical range of spin, and recent observations lend ev- 
idence that at least some Pop III stars had significant an- 
gular momentum. For instance, it has been found that some 
very metal-poor stars in the Milky Way (MW) halo have 
an a nomalous depletion of Li well below the Spite plateau 
(e.g. iFrebel et alj|2005l ; ICaffau et al. I l201ll ; iBonifacio et al.l 
12012? ). This is som etimes referred to as the 'meltdown' o f 
the Spite plateau (|Aoki et al.ll2009l ; ISbordone et al.ll2010l ). 
Such Li destruction may have occurred through rotation- 
ally induced mixing (e.g. IPinsonneault et al.ll2002i ; see also 
discussion and references in lAsplund et alj|2006l ). 

Further observations by IChiappinT et al.l f|201lh also 
lend evidence for rapid rotation of previous generations of 
massive stars. They find anomalous enhancement of Ba, La, 
and Y within the globular cluster NGC 6522. These elements 
may have been produced through an enhanced s-process in 
rapidly rotating massive stars. Observations of N/O and 
C/O abundance ratios in metal-poor stars in the halo as 
well as da mped Lyman-q sys tems , presented in works such 
as those of lSpite et ail |2005l ) and lPettini et all (|2008l ) . can 
also be more easily explaine d by enrichment from rapidly 
rotat ing massive stars (e. g. IChiappini et al.l |2006l; iHirschil 
120071 , see also discussion in iMaeder fc MevnetH2012l ). 

Current observations of O- and B-type stars in our 
Galaxy reveal that massive stars can indeed be rapid ro- 
tators, with rotation rates as high as several tens of percent 
of break-up speed, up to over 300 km s" 1 (e.g. lHuang fc Giesl 
120081 ; IWolff et all 120081 ). Lower-metallicity massive stars 
such as those within the Magellanic Clouds have also been 
found to have faster a verage rotation rate s than stars of 
higher metallicity (e.g. iHunter et al.l [iooil ). The environ- 
ments in which these stars formed, however, differs fr om that 
of Pop III stars (see, e.g.. IZinnecker fc Yorkell2007l ). While 
Pop III stars form in DM-dominated minihaloes, massive 
stars today form within non-DM-dominated giant molecu- 
lar clouds. If a fraction of small Pop III stars were ejected 
from their minihaloes, however, their accretion would be 
disru pted and they wo uld remain low-mass and long-lived 
(e.g. ICreif et all 1201 ll ). If these Pop III stars were later 
reincorporated into the Milky Way or one of its satellites, 
we may eventually directly observe them and their rota- 
tion rates. Identifying these as Pop III furthermore re- 
quires that they do not accrete too la rge an amount of 
metals from the enriched I SM (see, e.g.. IFrebel et al1l2009l ; 
Ijohnson fc Khochfarl l201ll ). However, such Pop III stars 
have yet to be found. To examine the range of possible ro- 
tation rates for Pop III stars, and thus their potential for 
various spin-dependent evolutionary pathways, we therefore 
utilize numerical simulations, initialized in their proper cos- 
mological context. 

In the work of IStacv et al.l (|201ll ), we estimated the ro- 
tation rate of Pop III stars from a cosmological simulation 
that resolved scales as small as 50 AU, finding these stars can 
attain very rapid rotation, potentially as high as their break- 
up speed. This was deduced from the rotational and thermal 
structure of the star-forming disk on the smallest resolv- 
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Figure 1. Profile of rotational velocity within the most massive protostar of each minihalo, shown with the solid lines and measured 
from the centers of the protostars. (a): Minihalo 1. (b): Minihalo 2. (c): Minihalo 3. (d): Minihalo 4. Dashed lines show i^Kep f° r 
comparison. Vertical blue dash-dot lines denote the photospheric surface of the protostar, R p . Vertical blue dashed lines show location 
of the hydrostatic surface of the star, i?« . Dotted red lines show example profiles of solid-body rotation. Note the significant rotational 
support within large portions of the protostars. The secondary peaks beyond the protostellar edge in panels (b), (c), and (d) are due to 
smaller protostars that are in tightly bound orbits around the main protostars. 



able scales, though the behavior of the gas as it reached the 
stellar surface and accreted onto the sta r could not be fol - 
lowed directly. The more recent study bv lGreif et alj (|2012l ) 
employed simulations which reached significantly increased 
resolution, down to scales of 0.05 R©, in four different mini- 
haloes. They found high rates of mergers between protostars, 
such that about half of the secondary protostars within a 
minihalo merge with the primary protostar. This was the 
first simulation to follow fragmentation and merging with 
this level of resolution, while previous pioneering simula- 
tions i^ffjiompjara^ 

fe.g. lOmukai fc Nishill 19981 ; iRipamonti et alj2002f l . or three- 
dimensional but unable to follow the gas evolution beyond 
the formation of t he initial hydrostatic core. We analyze the 
iGreif et al.l (|2012l ) simulations to obtain improved estimates 
of Pop III rotation rates and internal structure. Because 
these calculations resolved scales within the surface of the 
protostar beginning from several different cosmological re- 
alizations, we can now study the effect of mergers on the 
rotational structure, and the variation of spin for protostars 



within different host minihaloes. In Section 2 we give an 
overview of the numerical methodology used in the simula- 
tions. In Section 3 we present our results, assessing the main 
caveats in Section 4. We conclude in Section 5. 

2 NUMERICAL METHODOLOGY 
2.1 Initial Setup 

Our w ork utilizes the simulation output of IGreif et al.l 
(|2012l ), described in detail therein. Briefly, the initial con- 
ditions for these simulations were taken from IGreif et al.l 
|201lf ). in which the calcu lations were p erformed with the 
moving-mesh code AREPO (|Springe]||20roh . The original cos- 
mological DM simulations employed boxes of lengths 250 
and 500 kpc (comoving). They were initialized with 128 3 
and 256 3 particles at z = 99, assuming a A cold dark matter 
(ACDM) cosmology with fi m = 0.27, fi b = 0.04 6, h = 0.71, 
and u 8 ranging from 0.81 to 1.3. See table 1 of IGreif et al.l 
|2012f ) for details. The simulations are then followed until 
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Figure 2. Profile of radial velocity within the most massive protostar of each minihalo (solid lines), measured from the protostcllar 
centers, (a): Minihalo 1. (b): Minihalo 2. (c): Minihalo 3. (d): Minihalo 4. Vertical lines employ the same convention as in Figure[T] Note 
that « ra< ] is nearly zero within il*, indicating the hydrostatic equilibrium within the newly formed protostars. Dashed lines show M tnl b, 
which reaches supersonic levels outside of R p , but is negligible within the protostars. 

the formation of the first minihalo with virial mass greater 2000 AU is again extracted for further refinement, such that 
than 5 x 10 5 Mq. the final simulations resolve scales as small 0.05 R©. 



2.2 Cut-Out and Refinement 

Once the site of the first minihalo has been determined, new 
initial conditions at z = 99 are generated. Particles that 
will become part of the first minihalo, as determined by the 
original simulations, are replaced by 64 DM particles and 
64 mesh-generating points to be used in the hydrodynamic 
calculations. Cells and DM particles with increasingly large 
distances from the high-resolution region are replaced with 
cells and particles of progressively higher mass. Therefore, 
regions farther from the target minihalo are more coarsely 
resolved, reducing the total number of cells and particles to 
be followed for these refined simulations. 

After a cell in the refined simulations reaches a density 
of fO 9 cm -3 , the central 1 pc is extracted and used as the 
initial conditions for a further-refined simulation with reflec- 
tive boundary conditions. These new simulations are then 
followed to a density of fO 19 cm -3 , after which the central 



2.3 Chemistry, Heating, and Cooling 

The chemical and therm al network, described in detail in 
iGreif et alj (|201ll . I2012T ), follows the evolution of H, H+, 
H~, H+, H 2 , He, He+, H e ++ , D, D+, HD, and free electrons 
(|Glover fc Jappsenll2007l . Clark et al 2011a). Particularly im- 
portant chemothermal processes include cooling through H2 
rovibrational transitions at low densities, three-body H2 for- 
mation at densities above 10 8 cm -3 , and collision-induced 
emission (CIE) at densities greater than 10 14 cm -3 . H2 
collisional dissociation cooling also provides an off-set to 
compressional heati ng at high densities (e.g.. IOmukai|[2C}00l ; 
lYoshida et al]l2008l ). 

2.4 Extraction and Identification of Protostars 

Each protostar, along with its mass M* and photospheric 
radius R p , was identified using the procedure described in 
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d [R 8 ] d [R e ] 

Figure 3. Profile of the enclosed mass M cnc in the protostcllar regions of each minihalo, as measured from the centers of the protostars 
(solid lines). Also shown is the profile of f! in units of km s — 1 Rq (dashed lines). Vertical lines have the same meaning as in Figure [T] 
Note the flattening of M enc beyond the edge of each protostar, and the roughly solid-body rotational profiles inside the protostars. (a): 
Minihalo 1. (b): Minihalo 2. (c): Minihalo 3. (d): Minihalo 4. 



iGreif et al l (|2012l ). In short, each new protostar was found 
by searching for cells with density surpassing n — 10 19 cm" 3 . 
The densest of these cells was used as the center of the can- 
didate protostar. If the protostellar center was located out- 
side the radius of any previously identified protostars, it was 
counted as a new protostar. The protostellar boundary was 
defined as the edge of the photosphere, which was calcu- 
lated by determining the spherically averaged radius where 
the optical depth reached unity. 

Once a candidate cell was located, the optical depth 
At was determined within ~2x 10 6 bins. Bins were found 
by dividing the area into -/V ang ~ 10 4 uniformly-spaced an- 
gular regions that were further divided into JV rad = 200 
logarithmically-spaced radial segments between 0.01 and 10 
AU, centered on the candidate cell: 



(1) 



tegrated optical depth for each bin was then summed along 
each radius: 
i=k 

IGreif et alj |2012l ) next calculated the radial index fc cr i t at 
which the spherically averaged escape fraction, 



/V, 



1 - exp (-Tf,k) 



Tj,k 



(3) 



drops to /Scrit = 1 — exp (— 1) ~ 0.63, which corresponds to 
an optical depth of unity. The photospheric radius R p was 
then set equal to rk(k C Ht), while M» was given by the mass 
enclosed within R p . 

In our current study, we furthermore define a 'hydro- 
static' radius R„, defined as the distance from the protostel- 



where j denotes the angular region, k the radial segment, 
Pj t k the mass-weighted density of the bin, Kj t k the Rosse- 
land mean opacity, and Art the radial extent of the bin. 
Opacities were taken from iMaver fc Duschl (|2005l ). The in- 



value of ft. As will be seen in Section 3.2, the radius i?* < R p 
where k reaches a peak corresponds more precisely to the 
location of the protostellar accretion shock. 

We focus our study on four protostars, the most mas- 
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Figure 4. Radial profile of various properties of the most massive protostar in each halo, (a): Number density, (b): temperature, (c): 
opacity k and (d): radiative luminosity L T with respect to distance from the protostellar center. Black solid line represents the protostar 
from Minihalo 1, blue dotted line from Minihalo 2, green dashed line from Minihalo 3, and red dash-dotted line from Minihalo 4. Short 
vertical lines of corresponding style in panels (a) and (b) show the location of the photospheric surface R p . In panels (c) and (d), the 
vertical lines instead denote i?*. Thick yellow dotted lines represent example analytic stellar structure solutions. From bottom to top, 
these solutions are for polytropic indices of n =1.5, 3.0, and 4.0, with normalizations based on the mass and radius of the Minihalo 1 
protostar. Comparing the solid black line with the yellow dotted lines, the protostellar density is best described by the n = 3 solution, 
while the temperature is better-described by the n = 1.5 solution. In panel (d), upper thin lines correspond to L a g . As expected, L T 
approaches L a g near the protostellar surface. Thin-lined sections of the L r profiles represent regions where dT^/dr is positive, leading to 
negative L r values. In these sections we thus show the modulus of L T . Note that at these times the conditions for nuclear ignition have 
not yet been satisfied. 



sive proto star found within each of the four minihalos sim- 
ulated in iGreif et alJ (|2012h . This allows us to probe the 
early stages of rapidly accreting protostars within a range 
of host minihalo environments. Though there are a number 
of smaller protostars within each system, we here focus on 
the protostar from a given minihalo that is most likely to 
later become a massive star with the greatest impact on its 
surroundings. 



3 RESULTS 

3.1 Protostellar Rotation 

The significant rotational support within the protostars at 
the end of each simulation can be seen in Figure [1] where 



we show the radial profile of v rot . For each radial bin we 
determine its rotational velocity v rot as follows: 



Wrot = 



,(4) 



where rrii is the mass of a single hydrodynamic cell in the 
radial bin, v I0 t,sc,i is the x component of cell's rotational 
velocity, and M is the total mass within the bin. 

Figure [T] furthermore displays the Keplerian velocity 
t>Kc P = (GMenc/r) 1 , where M enc is the mass enclosed 
within distance r from the center of the protostar. Also 
shown is an example solid-body rotational velocity profile, 



V solid 



= rimaxr, where Omax = Wr, 



x/Rp, and t) Iot , max is 
the maximum value of v Iot found within R p . At the proto- 
stellar surface, the rotational support varies from ~ 80% of 
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Figure 5. Physical and adiabatic temperature gradients with respect to distance from the center of the protostar, shown at the final 
snapshot time tj} n of each simulation. Solid lines denote V p hyai an d dotted lines denote V at j. Dashed lines show Vn out to _R», while 
beyond this radius our approximation for Vn no longer applies since the gas is not undergoing solid-body rotation in these regions. 
Vertical dash-dotted lines show the location of the photospheric edge of the protostar at the end of the simulations, while vertical dashed 
lines denote the location of R*. Convection occurs where Vphy S > V a( j + VnsinS. Thus, these protostars are generally non-convective at 
these early times, where the conditions for nuclear ignition have not yet been satisfied. 



VKep for Minihaloes 1 and 3 to > 95% for Minihaloes 2 and 
4. In panels b, c, and d, tightly bound secondary protostars 
are evident as extra peaks in the v Iot profiles beyond the 
edges of the main protostars. 

For comparison, the profile of radial velocity u ra d can 
be seen in Figure [2] The large v la d outside R, signifies the 
radial inflow onto the protostellar surface, while u rac j falls 
to zero inside i?* due to the hydrostatic equilibrium of the 
gas. We additionally show the turbulent Mach nu mber M tur b 
over a range of radial bins in Figure [2] (see also iGreif et al.l 
l2012h . defined as 

i 

where c s is the sound speed of the radial bin, m% is the mass 
of a cell with index i contributing to the bin, and M is the 
total mass of the bin. Turbulence increases where there is 
radial inflow, and thus is approximately sonic outside of R, 
while becoming supersonic outside of R p . However, it is neg- 



ligible inside of R, where the gas has attained hydrostatic 
equilibrium. 

Within the protostar, vkc P roughly follows i>kc P oc r, 
the same as u so iid, due to the nearly constant-density pro- 
tostellar cores which yield M enc oc r 3 . The decline of VKep 
at radii beyond the protostar coincides with the leveling off 
of the enclosed mass M enc (Fig- EJ , at which point wk cp will 
scale approximately as r~ 1//2 (see also Clark et al. 20f fb). 

Along with «Kc P , rotation rates within the protostars 
also follow an approximately solid-body profile. This is ap- 
parent in Figure |3j which along with M Gnc also shows the ra- 
dial profile of the angular velocity £7 = v Iot /r, in units of km 
s _1 Rq 1 . The protostars generally exhibit a near-constant 
f2 ranging between ~ 5 x 1CF 7 to > 1.4 x 1CT 6 s _1 . Near the 
protostellar surface and beyond, f2 shows the expected drop- 
off as roughly f2 oc r -3 / 2 . Such a profile profile results from 
Wrot remaining at significant fractions of vkc P throughout 
the protostar, because gas from the rotationally-supported 
disk maintains much of its rotational support as it accretes 
onto the protostar. Note, however, that the S7 profiles are 
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Figure 6. Evolution of various properties of the most massive protostar in each halo, (a): Mass, (b): photospheric radius R p (thick lines) 
and hydrostatic radius ij* (thin lines), over time. Line styles have the same meaning as in Figure l4l 



not as flat as the inner p rofiles of models presented in, e.g., 
iMevnet fc Maed"etl (|2000t) . Deviation from a solid-body pro- 
file is due to angular momentum redistribution provided by 
torques from infalling gas and secondary protostars, allow- 
ing for some differential rotation such that Q — Q(r). 

For stars that are not undergoing accretion or mass loss, 
redistribution of angular momentum can also occur through 
processes such a s shear stress, meridiona l circulation, and 
convection (e.g.. IMevnet fc Maed"eill200Ch . These processes 
require timescales longer than the times followed in the sim- 
ulation. Eddington-Sweet circulations, for instance, occur in 
radiative regions of a star on the timescale tss- In a uni- 
formly rotating star this is approximately 



tES = 



GMl GAL 



GM, 

n 2 m 



(e.g., iMaeder fc MevnetJ [2012] ), where £kh is the Kelvin- 
Helmholtz timescale. M*, i?*, and L, are the mass, ra- 
dius, and luminosity of the protostar. Typical values for M* 
and R* at the end of this simulation are ~ 0.5 Mq and 
50 R,0, while R p is a few times larger (> 100 Rq). The 
mass-weighted value of the radiative luminosity L r within 
R p , for regions where L r is positive, ranges from ~ 10 _,i 
Lq to 1 L (see also following section). Average L r values 
within _R* are smaller and range from ~ 10~ 5 Lq to 10~' ? 
L0 since the simulations are still in early pre-nuclear igni- 
tion stages. If we use larger photospheric values to estimate 
Ies, 1 L-0 and 200 R,0, we then find typical tEs values of ~ 
1000 yr, well in excess of the simulated time. For the much 
smaller luminosities encountered in the interior even longer 
times would result. 

Note that for stars that have contracted to the main- 
sequence (MS), tES will typically be several orders of magni- 
tude longer than our estimate, and low-mass stars will have 
tes values even longer than their MS lifetime. Also note that 
our estimate for L* is several orders of magnitude lower than 
what might be expected in the outer regions of the star (see 
following section), or from the accretion luminosity L acc : 



GM*M 
~R* 



For M ~ 5 x 10~ 2 Mq yr -1 , our typical protostar will have 
iacc ~ 4000 Lq. These high luminosity values, however, are 
not representative of the radiative flux through the majority 
of the protostellar interior, so we use the lower values of L« 

quoted above. 

These simulations confirm the prediction of lStacv et alj 

(2011), based on approximate modeling of the sub-grid 
physics, that the significant rotational support found in the 
large-scale ~1000 AU star-forming disk down to a few tens 
of AU would also persist on unresolved AU and stellar scales. 
Furthermore, our analysis demonstrates that the high rate 
of protostellar merging does not lead to significant transfer 
of angular momentum away from the main protostars. 



3.2 Internal Structure 

In Figure 3J we show radial profiles of additional quantities 
for the most massive protostar in each minihalo, taken at 
the final simulation output. We compare the density and 
temperature structure of the protostar of Minihalo 1 to a 
range of polytropic models (yellow dotted lines in Fig. 
where 



p(r) = p c 9 n {r), 
and 

T{r) = T c 8(r). 



(8) 



(0) 



6(r) is determined through solving the Lane-Emden equa- 
tion, which describes the relation between 6{r) and the di- 
mensionless radius e(r) — r/r n depending u pon the poly- 
tropic index n (see, e.g., lHansen et al. 1 12004 ) The normal- 
ization factors r„, p c , and T c are determined by the stellar 
mass and photospheric radius, M* and R p : 



Rp/ei 



(I) B; 



(7) 



(10) 
(11) 

(12) 
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Figure 7. Rotational velocity w ro t,* over time for each of the main protostars. Dashed lines show the evolution of «Kep for comparison. 
Bottom lines of each panel also show the number of mergers undergone by the main protostar over time. The stars maintain a rotational 
velocity that is a significant fraction of ^Kcp> and the Minihalo 1 star even temporarily has v TO t > «Kcpi though this is an artifact of our 
spherical averaging scheme and the protostar's distortion during its initial merger event. 



(n + l) 



-eit 



/ G/im H \ f M» 
V k B )\Rm 



(13) 



where ei and 9^ are the values of e and d6/de at = 0. 

The mass of a hydrogen atom is mH, fcrs is Boltzmann's con- 
stant, and fi the mean molecular weight. We take the mean 
molecular weight to be n = 1.22, which is the average molec- 
ular weight of all cells with n > 5 x 10 19 cm -3 . We note that 
[i is not constant throughout the protostar, contributing to 
imperfect fits between the simulation and the polytropic so- 
lutions. However, the protostar of Minihalo 1 seems to best 
follow models with n ranging from 1.5 to 3. This applies to 
the main protostars of the other minihaloes as well, though 
we do not display these fits in Figure [3] 

The Rosseland mean opacity k , shown in panel fc) of 
Fig- HI is determined as described in lGreif et al.l |2012t ). The 
peak in each k curve coincides with the protostellar accre- 
tion shock, roughly corresponding to where the temperature 
falls below ~ 2 x 10 4 K and density below ~ 10 18 cm -3 , with 
the opacity rapidly dropping at radii beyond the shock. The 
re peaks mark where the predominant opacity contribution 
transitions from H bound-free (bf) absorption to H~ bf ab- 



sorp tion as the gas phase c onverts from ionized to neutral 
(e.g. iMaver fc Duschj 120051 ). At larger radii beyond ~ 100 
R,0, the main contribution to opacity is from H2. 

Panel (d) of Fig. [4] estimates the radiative luminosity 
profile L r of each star, defined as 



L r (r) 



47ra cr 2 dT A 
3np dr 



(14) 



where a is the radiation constant, r is the distance from the 
stellar center, and dT^/dr is taken directly from the simula- 
tion output. The above equation describes the luminosity at 
a given radius using a diffusion approximation to radiative 
transfer and assuming that all energy is transported through 
the protostar by radiation. In this case the protostellar lu- 
minosity L t is equal to L r . The simulation, however, did 
not include full modeling of radiative transfer, and instead 
assumed a local escape fraction using the Sobolev approxi- 
mation. The resulting L r profile should be interpreted only 
as an approximation. For comparison purposes, panel (d) of 
Figure U also shows L c g(r) = 47rr 2 <TSBT 4 ('~). As expected, 
for each protostar L r approaches L e g(r) near r = R», where 
T approaches the effective temperature T e g. We also point 
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Figure 8. Level of rotational support /3 over time for each of the main protostars (solid lines). Dotted lines show the level of rotational 
support of the gas surrounding the protostar out to distances of 10i? p . Dashed lines show the ratio of v TO t to 'UKep at the protostellar 
surface R p . After the first ~ 1 — 2yr, /3 fluctuates around a nearly constant value of 0.3. 



out that regions of the star where dT 4 /dr is positive lead 
to negative values of L r . Regions with negative L r values 
correspond to the thin-lined sections of the profiles, where 
instead of the negative L r values we have plotted the mod- 
ulus. 



The accretion flow outside of the photospheric surface 
R p is subject to supersonic turbulence (see Fig. [2]), which 
provides an additional source of heating. The values of L r 
just beyond R p are typically on the order of 1000 Lq , similar 
to the value L acc ~ 4000 Lq quoted in the previous section. 
The main contribution to the luminosity is thus the impact 
of accreting gas onto the protostellar surface, though there 
is some contribution from the turbulence generated by the 
inflow as well. For M tur b = v t urb/c s ~ 2, and c s ~ 5 km 
s _1 , the turbulent velocity v t urb is 10 km s _1 . Referring to 
Figure |31 the turbulent region beyond R p (e.g., the region 
between 200 and 300 Rq for the protostar of Minihalo 2) 
contains a small mass of approximately 0.01 Mq, yielding a 
kinetic energy of 10 43 erg. With typical radial velocities of 3 
km s _1 , the timescale for gas to cross the turbulent region 
is ~ 10 7 s, yielding an energy production rate of 10 36 erg 



s _1 , or ~ 250 Lq. While not negligible, this is significantly 
smaller than L acc 

We furthermore point out that the effects of radiation 
pressure were not accounted for in the simulations. We ex- 
pect radiation pressure from direct ionization to be a negligi- 
ble effect, as the radially distended protostars have a typical 
effective temperature of ~ 3000-4000 K, assuming they have 
surface luminosities of L* = L acc ~ 4000 Lq. Note that this 
agrees well with the gas temperature at the photospheric 
surface of the protostars, ~ 3000 - 5000 K (see panel b of 
Fig. |4]). These protostars thus do not yet ionize their sur- 
rounding gas. 

For gas within the protostar, we can estimate the radi- 
ation pressure provided due to L r using 



K,pL r 



dr 4:Tvr 2 c 
which can be further simplified to 



Prad ~KpAr 



47rr 2 c 



47rr 2 c 



(15) 



(16) 



Radiation emitted from just within the photospheric surface, 
usually at r ~ 100 Rq , has typical L r values of 1 Lq , and by 
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definition r ~ 1. This yields P ra d ~ 10" dyne cm" . The 
radiation pressure exerted upon the outer atmosphere of the 
protostars is thus negligible compared to the thermal pres- 
sure, Ptherm = nksT ~ 5 x 10 4 dyne cm -2 for representative 
values of n = 10 17 cm" 3 and T = 4000 K. 

Inside the protostar, the average L r within R* has lower 
values of ~ 10" 4 Lq. However, r is orders of magnitude 
larger, r ~ npR, ~ 10 11 for k = 1500 cm 2 g _1 , R, = 50 R Q , 
and p — 3 x 10 -5 g cm -3 . With these values we can approx- 
imate that radiation within the protostellar accretion shock 
exerts a pressure P ra( j ~ 10 4 dyne cm -2 . On the other hand, 
typical gas temperatures and densities in the inner protostel- 
lar regions are T ~ 3 x 10 4 K and n ~ 3 x 10 19 cm" 3 , yielding 
a thermal pressure of P t hcrm ~ 10 8 dyne cm -2 , again signif- 
icantly greater than the radiation pressure. 

We can also estimate the strength of radiation pressure 
due to Thomson scattering by comparing L r within the ion- 
ized region of the protostar to the Eddington luminosity, 
£<Edd = 47rGM*mHc/(TT, where or is the Thomson scatter- 
ing cross section. Our typical protostar has Z/Edd ~ 10 4 L0, 
orders of magnitude higher than L r within the ionized gas. 
Radiation pressure will thus not be a significant effect until 
nuclear burning has commenced. 



3.3 Convective Instability 

In Figure[S]we estimate the regions of the protostar that are 
unstable to convection, determined by the relation between 
Vphys and Vad- As stated by the Schwarzschild criterion, 
convection in non-rotating stars is expected if V p h ys > V a d, 
where 

Vphys =(^-5) (17) 



, d In P ) phys 

is the actual temperature gradient within the protostar 
which we estimate directly from the simulation output. Fur- 
thermore, 



/dlnT\ 
\d\up) t 



(18) 



is the adiabatic variation of temperature with pressure. 

To account for rotation and the restoring effect of 
angular momentum, we must instead use the related 
Solberg-Hoiland criterion fo r convective instability (e.g. 
iKippenhahn fc Weigertlll990h : 



V phys 

Here 



> V a 



Vnsini? . 



(19) 



(20) 



v , _ H p 1 d(i]V) ^ 

where Hp is the pressure scale height, (/ grav the gravitational 
acceleration, $ the colatitude angle from the axis of rotation, 
u) — rsim? the distance to the rotation axis, and 



5 = - 



dlnp 



(21) 



.dhiT/p 

For solid-body rotation which very roughly describes our 
protostars within R t , this expression for Vn can be simpli- 
fied to 

n 2 H a 



Vn =4 



<7gr a v 8 



(22) 



We approximate t hat H v ~ -P/(fl E rav P), where P is the gas 
pressure (see, e.g-. lMaeder et al ] |2008h . For simplicity, m our 
estimates we also set sim? equal to one. 

If all energy transport is through radiation, we may also 
write L, — L r and Vphys = V ra d, where 



Vr E 



/dlnT\ 
\dlnPJr 



Pk u 



ad WnacG T 4 M e; 



(23) 



To account for the effects of rotation, we may furthermore 
replace M enc with 



M r = Me, 



1 - 



(24) 



von Zeipellll924h where po is the average stellar density and 
is typically ~ 10" 5 g cm" 3 within _R*. For our protostars the 
parenthetical factor in the above equation is of order unity 
and ranges from 0.8 to 0.9. 

Figure [5] shows V p h ys , V a d, and Vn at the end of 
each simulation. Note that the estimate of Vn begins to 
break down outside of 7?*, where Q starts to deviate from 
a solid-body profile, so we only show Vn out to this radius. 
Throughout most of the protostar, V a d is the largest term 
and prevents convection within the protostars. Vphys and 
Vn instead maintain values closer to zero with the excep- 
tion of regions within ~ 10 -Rq. We emphasize, however, 
that improved simulations, including the proper treatment 
of diffusive transport of radiation, will be necessary to ac- 
curately model convection within Pop III protostars. Bear- 
ing this caveat in mind, at the end of the simulations only 
one protostar, that of Minihalo 2, exhibits possible central 
convection. Though not shown in Figure [5j at earlier times 
other protostars also have short periods where the convec- 
tion criterion is satisfied, giving some indication that Pop 
III protostars may have typically undergone periods of con- 
vection in their cores as they grew in mass. However, overall 
the protostars appear non-convective. 

Prototars on the Hayashi track are expected to have 
convective zones, but if the protostar is growing then 
the presence of an accretion shock will return it to sta- 
bility against convection (e.g. Starrier 1988b, see als o 
IWuchterl fc Klessenl l200ll ; IWuchterl fc Tscharnuter] 120031 ). 
However, previous one-dimensional m odeling of Pop HI 
evolution under rapid accretion (e.g. 



lOmukai fc Pallal 120031 ; iHosokawa et al.l 



Stahler et ail Il986l ; 
20ld ) indicate that 



growing protostars can still become convectively unsta- 
ble after the o nset of nuclear burning. In particular, 
IHosokawa et al.l (|2010h find that metal-free protostars of 
mass < 1 Mq accreting through a disk may briefly have 
a central convective zone after deuterium burning begins. 
However, this zone soon disappears while a longer-lasting 
surface convection zone emerges as the protostar grows and 
the D-burning region moves to the outer layers. For spher- 
ical accretion, however, deuterium burning and convective 
zones do not appea r until the protostars h ave attained higher 
masses (< 10 M^. lQmukai fc Pallal l2003h . 

Without modeling of this nuclear ignition or of radiative 
transfer, however, we do not expect sustained convection to 
appear in our simulations. We furthermore note that, despite 
the differences in opacity between Pop III, low-metallicity, 
and solar-metallicity protostars, appearance of convection 
zones associated with the onset of deuterium burning is ex- 
pected regardless of metallicity. In all cases the main contri- 
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t = 8.31 yr I t = 11.38 yr 




1 AU (215 R s ) I 1 AU (215 R«) 



Figure 9. Morphology of the main protostar within each minihalo, where images arc oriented to show the face of the disk. Each panel 
has a width of 5 AU. A line of length 1 AU is also shown in each panel for comparison. Dashed white circles depict the size of the 
photosphere of the main protostars, while solid white circles denote R*. The color scale ranges from density of 10 12 cm -3 to 10 21 cm~ 3 . 
The diskiness and spiral structure of the accreting gas is readily apparent. 



butions to opacity are not metals bu t H absorption an d 
photoionization (e.g. Stahler 1988a; iDurisen et alj Il989h . 
though the extra contribution to the electron fraction from 
metals d oes still lead to higher opacities in solar-metallicity 
gas (e.g. lAlexander et al 

3.4 Time Evolution 

Figures [6] and [7] show how various properties of the proto- 
stars evolve over time. Mass M* and radii R D and R, of 
the p rotostars are determined as described in iGreif et al.l 
(|2012h and Section 2.4. In Figure [7] we determine u ro t,* 
similarly to the method described in Equation 4, except 
that the sum is performed over all hydrodynamic elements 
within R p , and M is replaced with M*. When compared 
with vkc P ,* = (GMt/Rp) 1 ^ 2 , we see that each of the stars 
maintains a rotational velocity that is a significant fraction 
of its break-up speed. The protostar of Minihalo 1 (panel a 
of Figure [7| even appears super-Keplerian for a brief period 



around 4 yr. However, this is an artifact of our spherical 
averaging scheme, and the super-Keplerian period coincides 
with when the protostar is undergoing its first merger and 
exhibits a shape that is more flattened and bar-like than 
spherical. We here briefly mention that many of the sec- 
ondary protostars in each minihalo also maintain very large 
rotation rates. The majority of secondary protostars have 
values of v r ot,*/i>Kep,* ranging from ~ 70-90% at the final 
simulation outputs, though some have lower values of ~ 20- 
30%. 

We also quantify the level of rotational support with 
the parameter /3, which we define as the ratio of the specific 
angular momentum of the protostar, j„, to its Keplerian 
angular momentum jKep ■ We thus have 

P=j*/3Kep, (25) 

where 
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11.38 yr 




1 AU (215 R.) 



Figure 10. Same as previous figure, except with images oriented to show the disk edge-on. Each panel has a width of 5 AU. Note the 
flattened structure of the protostars at these early times. 



J. = 



1/2 



(26) 



The index i refers to values for individual gas elements 
within the protostar, and 

jKcp = (GM t R p ) 1/2 . (27) 

In Figure[8l it can be seen that each of the stars remains sub- 
stantially rotationally supported and maintains a relatively 
constant value of /3 ~ 30%, with some fluctuation above and 
below this value. The variations in /3 also tend to follow the 
fluctuating rotational support in the surrounding disk, out 
to 10-Rp (dotted line in Figure [5]), with the exception of the 
Minihalo 1 protostar. In addition, the ft oscillations roughly 
follow the fluctuations in the ratio of rotational to Keplerian 
velocity at the protostellar surface (dashed line in Figure [SJ. 
Each protostar generally maintains v rot of > 80% of «Ke P at 
its surface, demonstrating that the rapid protostellar rota- 
tion is set by the nearly Keplerian rotational speeds at the 
boundary between the inner disk and protostellar surface. 
The rapid rotation at the inner disk boundary in turn orig- 



inates from the angular momentum of the larger-scale disk. 
This Kep lerian inner disk is likely to persist to later times 
(see, e.g. JStacv et al.ll201lh . Though these simulations only 
follow the initial stages of protostellar accretion, as the pro- 
tostars grow by an order of magnitude their high spins may 
be maintained through continued accretion from a nearly 
Keplerian disk. 

This significant rotation causes the protostars to have 
flattened structures, as is visible in Figures [9l and [TOl where 
we depict the density structure around the main protostars 
at the end of each simulation. Note that the flattening is 
much more significant on scales of R p than on the smaller 
scales of i?* . Furthermore, the main protostars of Minihalo 
2 and 4 exhibit even more asymmetry within R p than those 
of Minihalo 1 and 3. This is because in the former case, the 
protostars are undergoing mergers at the t imes shown, lead - 
ing to greater flattening. As discussed in iLin et alj (|201lf ). 
such deviations from spherical symmetry in stars undergoing 
disk accretion are to be expected as angular momentum is 
exchanged between the star and disk through gravitational 
torques. Similar to the stellar rotational evolution seen in 
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their calculations, each protostar rapidly spins up over the 
first few Keplerian rota tional periods tk cp , where in our case 
typically tk cp ~ 1 yr. iLin et all (|201ll ) furthermore find a 
longer-term stabilization of the spin evolution, where the 
stellar spin may remain at ~ 50% of the break-up value, or 
may undergo a slow decline, depending on the properties of 
the surrounding disk. Their simulations address only non- 
fragmenting disks. Nevertheless, even given the fragmenta- 
tion and protostellar merging in our simulations, we still see 
a similar lack of significant evolution in /3 for any of the 
protostars. 



3.5 Dependence on Minihalo Characterstics 

The initial collapse of the minihalo gas f rom IGM to proto- 
stellar densities was discussed in detail in lGreif et alj (|2012T ) . 
They found that the radial and rotational velocity structure 
of each minihalo, as well as density and temper ature evolu- 
tion, were consistent with previous studies (e.g . Abel et all 
120021 : iBromm et all |2002l : lYoshida et ail 120061 l2008h ■ Fur- 
thermore, the gas collapse within each minihalo generally 
follows the Lars on-Penston so l ution for isoth ermal and self- 
gravitating gas (|Larsodfl96i iPenstonl fl969h ■ The angular 
momentum profiles of the central few hundred solar masses 
of each minihalo are all the same to within a factor of 
a few, in addition to showing similar ag reement to previ- 
ous cosmological simulations (Figure 1111 lAbel et al.l 120021 : 
lYoshida et al]|2006l : Istacv et alj|20ich . 

The minihaloes are not perfectly identical, however, and 
it is interesting to see whether variation in properties of 
the host minihaloes lead to any systematic variations in the 
properties of their protostars. Figure [T2"l shows the final w ro t,* 
of each star with respect to the spin parameter, A, and mass 
of the minihalo. Also denoted is the typical spin parameter 
of DM haloes, A = 0.05, as measured in large-scale cos- 
mological TV-body simulations ()Barnes fc Efstathioul 1 19871 : 
IJang-Condell fc HernqrustboOll l. where A is defined as 



J\E 



1/2 



GM, 



5/2 



(28) 



Here J, E, and Mhaio are the total angular momentum, en- 
ergy, and mass of the halo, respectively. We have too few 
sample minihaloes, however, to find an obvious correlation. 
Furthermore, when the protostars have first formed, as well 
as at the end of the simulations, the rotation axes of the 
protostars are not aligned with that of their minihaloes. The 
angle between the rotation axes varies from 30° to 160° . This 
further demonstrates a lack of connection between the spins 
of the protostars and their host minihaloes. We additionally 
checked for relationships between the total number of proto- 
stellar fragments formed, Nf rag , and A and w ro t,*- Again, we 
did not simulate a sufficient number of minihaloes to find 
statistically significant correlations, though such an inves- 
tigation with an increased number of minihaloes would be 
worthwhile for future work. 



4 CAVEATS 

We point out the caveat that magnetic fields, which were 
not included in these computations, may affect the angular 
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Figure 11. Angular momentum profile of gas within each mini- 
halo. Each minihalo is denoted by the same line style as in previ- 
ous figures. Thin yellow lines denote angular momentum profiles 
found in separate cosmological simulations. Solid yellow line is 
taken from Stacy et al. (2010), dotted yellow line from Yoshida et 
al. (2006), and dashed yellow line from Abel et al. (2002). Thick 
black dashed line shows an approximate powerlaw fit to these 
profiles, j oc M cnc . The profiles are all very similar even for a 
variety of cosmological realizations. 



momentum build-up of the protoste llar disk and the subse- 
quent protoste l lar ro t ation rate (e.g . iTan fc Blackmar]|2004l; 
Silk fc Langerl l200rj; iMaki fc Susal l2007t ^Schleicher et al l 



20101 ; iFederrath et alj|2011l : ISur et alj|2012l ; Schober et al. 
2012a, 2012b). In addition, magnetic fields may 'disk-lock' 
the star to a certain rotation rate that depends upon its 
mass , accreti on rate, magnetic field, and radius (iKoenigll 
Il99ll ; see also lShu et al.lll994 iMatt fc Pudritj|2005h . 

Magnetic fields may also allow for the operation of the 
Spruit-Taylor dynamo within the protostar (|SpruirJl2002l ') . If 
this mechanism is active in Pop III stars, it can facilitate 
angular momentum transport from the core to the enve- 
lope. This allows the star to maintain solid-body rotation, 
spinning up its outer layers and allowing the surface to reach 
critical rotation speed earlier in its evolution as compared to 
models which do not include magnetic torques. This helps 
to explain the difference in the results found by, e.g., Ek- 
strom et al. (2008a) and Yoon et al. (2012). For stars that 
undergo SN deaths, magnetic torques may additionally lead 
to a rotation rate of the collapsing iron core that is over 
an ord er of magnitude sm aller than those without magnetic 
fields (|Heger et al.l |2005l) . To model the formation of Pop 
III stars with further improved accuracy, future simulations 
shoul d include the effects of magnetic fields (e.g. iTurk et al.l 
120121) . 

Another important caveat to bear in mind is the short 
timescales followed in the simulation. The simulations fol- 
lowed ~ 10 yr, which is a small fraction of the Kelvin- 
Helmholtz time, £kh ~ 10° yr for massive stars, the typical 
timescale for stellar assembly. The rotational structure of a 
protostar only years after it has newly formed will not be 
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Figure 12. Comparison of the properties of the largest protostars with properties of their host haloes, (a): Initial and final measured 
rotational velocities of the largest prototar, i> r ot,*, versus the spin A of the host halo. Initial values arc in blue, while final values are in 
black. Red line shows the average A value for DM haloes (A = 0.05, e.g. Barnes & Efstathiou 1987). (b): Initial and final i> ro t,* versus 
minihalo mass M^aio- (c): Total number of fragments -/Vf rag formed with respect to spin of the minihalo. (d): Final « r ot,* versus total 
number of mergers with the main protostar N mcI g C . The diamond represents Minihalo 1, the triangle Minihalo 2, the square Minihalo 3, 
and the circle Minihalo 4. 



the same as that when the protostar enters or later leaves 
the MS. Longer simulations will be required to see how the 
stellar spin as well as the inner accretion disk will evolve. 
Only if mass infow continues to maintain high rotational 
support in the inner disk, and if the inner disk continues to 
transfer angular momentum onto the star, can large stellar 
spin rates be maintained. 



5 DISCUSSION AND CONCLUSION 

We use the cosmological simulations of iGreif et alj (|2012t ). 
which resolve four separate minihaloes down to sub-stellar 
scales (0.05 Rq), to analyze the rotation and internal struc- 
ture of Pop III protostars. This is the first simulation to pro- 
vide a direct view of Pop III protostellar structure from cos- 
mological initial conditions. We find the protostars quickly 
develop a roughly solid-body rotation profile, while their 
surface rotation velocities range from ~ 80 — 100 % of t* ep 
at the end of each simulation. Each of the four protostars 



examined maintains high rotation velocities even after un- 
dergoing multiple merger events. 

The protostars generally seem non-convective over the 
timescales simulated. However, future calculations which in- 
clude a proper prescription for radiative diffusion, as op- 
posed to the simplified escape probablity formalism em- 
ployed in Greif et al. (2012), will be necessary to properly 
model convection in Pop III protostars. There is also little 
evidence of correlation between the properties of each host 
minihalo and the spin of its largest protostar or the total 
number of protostars formed in the minihalo. More mini- 
haloes would be necessary to derive meaningful statistics, 
though it would be very informative for future simulations 
to more thoroughly examine correlations between minihalo 
properties and the Pop III systems they host. 

We furthermore note that future work will follow the 
protostellar evolution for significantly longer timescales. 
This will allow for a more direct determination of how the 
protostars evolve, and how this evolution is affected by ro- 
tation. As the protostars continue to grow and continue on 
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to the MS, rotation could alter the protostar's life in a num- 
ber of ways. Rotation rates which persist at sufficiently high 
speeds, for instance, may allow for mass loss through stellar 
winds generated at the 'fir" limit (see discussion in, e.g., 
iMaeder fc Mevnetl l20"il ). even though mass loss through 
line-d riven winds is expected to be minimal (e.g. iKudritzkl 
l2002h . This would reduce the final mass of the star and thus 
may alter the stellar death it will undergo. 

As described above, the metal production during the 
lifetime of the star wo uld also be gener ally enhanced (e.g. 
Ekstrom et al. 2008a. lYoon et alj|20l3 l The temperature 
and luminosity of the star will also be altered, and possi- 
bly greatly enhanced if r otational mixing is sufficient for 
the star to undergo CHE (|Yoon et alj|2012l ). though not all 
studies find that CHE can take place in rotating Pop III 
stars (e.g. Ekstrom et al. 2008a). CHE may furthermore 
provide a mechanism for a Pop III star to become a WR 
star and eventually a GRB without being in a tight binary 
dYoon fc Langeri 120051 ; IWooslev fc Hegeri l200rj ; lYoon et all 
120121 1 We finally note that CHE may also lower the min- 
imum mass at which a st ar will undergo a PISN death 
from 140 M fl to ~ 64 M (|Chatzopoulos fc Wheeler! |2012| ; 
lYoon et alj|2012f ) 

The high spins seen in these protostars for a range of 
minihalo characteristics already demonstrate that a Pop III 
star of a given mass will have a number of possible evolu- 
tionary pathways depending upon its rotation rate. Com- 
putational power is progressing to the point where three- 
dimensional simulations can begin to examine not only the 
mass growth but also the spin evolution of Pop III stars. 
Observational evidence will also provide constraints on the 
rotation rate of previous stellar generations. This evidence 
will most likely be the chemical signatures left in the old- 
est stars, though any putative low-mass Pop III stars that 
were ejected from their minihaloes may also someday be 
directly detected. Continually improving numerical and ob- 
servational efforts will allow us to probe the early Universe 
in ever-greater detail. 
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